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Abstract 

Introducing new physically motivated ansatze, we explore both analytically and nu- 
I merically the classical and absolute stabilities of a single Q-haW. in an arbitrary number 

' of spatial dimensions D, working in both the thin and thick wall limits. 

(N : 

; 1 Introduction 

^ ■ In a pioneering paper published in 1985 [1], Sidney Coleman showed that it was possible for a 

^ new class of non-topological solitons to exist within a self- interacting system by introducing 

. the notion of a Q-ball . His model had a continuous unbroken global U(l) charge Q (for reviews 
p • see [HI m O E]), which corresponds to an angular motion with angular velocity lo in the U(l) 
internal space. The conserved charge stabilises the Q-ball, unlike the case of topological solitons 
whose stability is ensured by the presence of conserved topological charges. Once formed, a Q- 
^ . ball is absolutely stable if four conditions are satisfied: (1) existence condition [T] - its potential 
^ . should grow less quickly than the quadratic mass term, and this can be realised through a 
. number of routes such as the inclusion of radiative or finite temperature corrections to a bare 
^ . mass, or non-linear terms in a polynomial potential [HE], (2) absolute stability condition - 
Xf^ . the energy Eq (or mass) of a Q-ball must be lower than the corresponding energy that the 
^ , collection of the hghtest possible scalar particle quanta (rest mass m) could have, (3) classical 
O ! stability condition [2] - the Q-ball should be stable to linear fluctuations; with the threshold 
j> ! of the stability being located at the saddle point of the D-dimensional Euclidean action, the 
^ ! bounce action 6*0;, [H], (4) fission condition pj - the energy of a single Q-ball must be less than 
the total energy of the smaller Q-balls that it could in principle fragment into. It turns out 
that for each of the four conditions to be satisfied we require: 

^ ^ oj dQ d'^S^ d'^En duj , , 

._<H<^.. S«<»Q. -^<0«^>0. ^<0«^<0 (1) 

where lo^ are the lower and upper limits of u; that the Q-ball can have. The lower limit, lo ~ a;_, 
can define thin wall Q-balls, either without [1] or with P, [10] the wall thickness being taken 
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into account, while the upper hmit, ~ 0;+, can define thick wall Q-balls in [TT] which may be 
approximated by a simple Gaussian ansatz [T^ . 

There is a vast literature on non-topological solitons, including Q-balls. They have been 
seen to be solutions in Abelian gauge theories [10], [IH [HI [15l [lEl [TT], in non-Abelian theories 
[m [la ED], in self-dual (Maxwell-) Chern-Simons theory [211 ESI ESI El] , in noncommutative 
complex scalar field theory [25], in models which include fermionic interactions [131 El IISI [13 
[2Sl [2Z], as well as in the presence of gravity [2HI [22] • Q-balls themselves have been quantized 
either by canonical [2] or by path integral schemes [SUl [SH [S2] • With thermal effects, it has 
been shown that Q-balls coupled to massless fermions are able to evaporate away [33] ; however 
at sufficiently low temperatures they become stable, and indeed they then tend to grow [HI [35]. 
The authors in [21 [lOl [361 [SZ] have discussed and analysed the spatially excited states of Q- 
balls, including radial modes as well as spatially dependent phase excitations. A more general 
mathematical argument concerning the stability of solitary waves can be found in [3HI [M]- A 
related class of objects to Q-balls are known as oscillons [iQl [lH [121 [131 [lH [15] or as I-balls [16], 
and recent attention has turned to the dynamics of these time dependent non-linear meta-stable 
configurations [HI [IHl [IH] • 

Standard Q-balls exist in an arbitrary number of space dimensions D and are able to avoid 
the restriction arising from Derrick's theorem [50] because they are time-dependent solutions. 
A few examples include polynomial models both for = 3 [HH [52] and for arbitrary D |12j : 
models with supersymmetry broken by gravity mediation [52]; and models with supersymmetry 
broken by gauge interaction [321 E31 [SI]- Returning to the case of D = 3, phenomenologically, 
it turns out that the Q-balls present in models with gravity mediated supersymmetry breaking 
are quasi-stable but long-lived, allowing in principle for these Q-balls to be the source of both 
the baryons as well as the lightest supersymmetric particle (LSP) dark matter particle [55]. On 
the other hand, Q-balls in models of gauge mediated supersymmetry breaking can be a dark 
matter candidate as they can be absolutely stable [H]. Both types of Q-balls have been shown 
to be able to provide the observed baryon-to-photon ratio [3^ . 

The dynamics and formation of Q-balls involve solving complicated non-linear systems, 
which generally require numerical simulations. The dynamics of two Q-balls in fiat Minkowski 
space-time depends on parameters, such as the relative phases between them, and the relative 
initial velocities [7l[S31[SS]- In addition, the main formation process through the Affleck-Dine 
mechanism [SZ] has been extensively examined in both gauge mediated [SB], gravity mediated 
[591 [nni [6l] , and running infiaton mass models [62]. As analysing individual Q-balls is difficult in 
its own right, it is extremely challenging to deal with multiple Q-balls. A number of analytical 
approaches to address that issue have been made over the past few years, e.g. [H31 [Ml [HS]- 
Multiple thermal Q-balls have been described in a statistical sense in [HUl [M]. 

In this paper we aim to analytically and numerically address stationary properties of a single 
Q-ball in an arbitrary number of spatial dimensions. The work will draw on earlier work of 
Correia and Schmidt [30] who derived analytic properties for the thin and thick wall limits of 
Q-balls in D = 3. Recently, Gleiser and Thorarinson [12] proved the absolute stability for thin 
wall Q-balls using the virial theorem. We will generalise the main results of [121 130] to the case 
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of arbitrary spatial dimensions, and in doing so will both analytically predict and numerically 
confirm the unique values of the angular velocity uJa for the absolute stability of the Q-balls via 
both the virial relation and thin wall Q-ball approximations. 

Moreover, we will obtain the classical stability conditions for the thin and thick wall ap- 
proximations, and discover the connections between the virial relation and thin or thick wall 
approximation for the characteristic slopes Eq/ujQ. In a companion paper [HTj we will inves- 
tigate dynamical properties of multiple Q-balls, including thermal effects, and their formation 
(see [68] for movies showing the dynamics of single Q-balls). 

This paper is divided into the following sections. In section [2] we review the properties of 
a Q-ball in an arbitrary number of spatial dimensions including the existence, and stability 
conditions. By introducing a number of different ansatze, in section O we present a detailed 
analysis of the solutions in the thin and thick wall limits. We then demonstrate the advantages 
of using two particular modified ansatze in section H] where we present detailed numerical results 
for the case of both degenerate and non-degenerate underlying potentials. Finally we conclude 
in section 



2 Q-ball in D dimensions 

We shall begin with a standard Q-ball ansatz which satisfies a Laplace equation called 
the Q-ball equation, and we will introduce the Legendre relations [2] which will make some 
computations easier. The existence of Q-balls places constraints on the allowed form of the 
potential, and introduces limiting values of cu, i.e. uj±, near which we may describe the Q- 
balls analytically using either a thick or thin wall approximation. We will then introduce three 
conditions for Q-balls to be stable [3]. Finally we will obtain the characteristic slope Eq/ujQ 
and minimum charge Qmin^ and propose approximate values for oja, the limiting frequency for 
absolute stability, using a virial theorem and showing that it does not rely on detailed analytic 
profiles and potential forms. 

2.1 Q-ball ansatz 

We consider a complex scalar field (j) in Minkowski spacetime of arbitrary spatial dimensions D 
with a U(l) potential bounded by f/(|0|) > for any values of 0: 

S = j d^^^x^ £, (2) 

where C = -^(^^'^S^^t^,^ - f/(|0|). (3) 

The metric is ds^ = g^^dx^dx^ = —dt^ + hijdx^dx^ and g is the determinant of g^y where 
yU,, V run from to and i, j denote spatial indices running from 1 to D. Now using the 
standard decomposition of in terms of two real fields = cre*^, the energy momentum tensor 
T^ij = + {symmetrising factors) and the conserved U(l) global current j^,u{i) via 
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the Noether theorem, we obtain 

T^u = {d,.ad,a + a^d,,edJ)+g^X, (4) 
ji,,u{i) = (r'^d^e. (5) 

Using a basis of vectors {f^l^a)} where nl'^j is time-hke and n^*^^ are space-hke unit vectors oriented 
along the spatial z-direction, the above currents give the definitions of energy density pE, charge 
density pq, momentum fiux Pi and pressure p: 

PE = T^un'"^^^n'(^y, pq = j^^n'^^y, P, = T^^nf^^nf^^; p = T^^n^^^nl^. (6) 

Defining the D dimensional volume Vd bounded by a (D — l)-sphere, the Noether charges 
(energy, momenta, and U(l) charge) become 



E= PE, Pi= P,, Q= Pq, (7) 

JVd JVd JVd 

where jy^ = J d^xy/h. Minimising an energy with a fixed charge Q for any degrees of freedom, 
we find the Q-ball (lowest) energy Eq by introducing a Lagrange multiplier uj and setting 
nr = (-l,0,0,...,0): 



Eq = E + uj[Q- pqJ, (8) 

= uQ + j^ (^^l^a^ + a\e-uy + {Vaf + a\Vey]+U^y (9) 

= ujQ + S^ (10) 

where U^j = U — ^cu'^a'^, a = ^ etc... and to will turn out to be the rotation frequency in the 
U(l) internal space. The presence of the positive definite terms in Eq. ([9]) suggests that the 
lowest energy solution is obtained by setting & = = 6 — uj = V6. The Euclidean action S^^ 
and the effective potential in Eqs. [TUl) are finally given by 

S.= [ ^(Va)2 + f/^, t/,^f/-Lv. (11) 
Jvd ^ ^ 

The second term in comes from the internal spin of the complex field. Following Friedberg 
et. al [2], it is useful to define the functional 



(12) 



where \'jJ^ is the corresponding Lagrange multiplier and / = J^^ . 

Given that the spherically symmetric profile is the minimum energy configuration [02], we 
are lead to the standard stationary Q-ball ansatz 

= a(r)e^'^*. (13) 



4 



Substituting Eq. (fT3l) into Eq. ([6]), we find 

PE = \(t'^ + f/ + pQ = Lja\ (14) 

Pr = V^-f/c., Pi = (15) 

where a' = ^ and is a radial pressure given in terms of tlie radially oriented unit vector 
= (0, 1, 0, . . . , 0). Without loss of generality, we set both u and Q as positive. 



2.2 Legendre relations 



It is sometimes difficult to compute Eg directly, but using Legendre relations often helps [2|[HU]. 
In our case, from Eq. ffTUj) and Eq. ffT^ we find 



dE, 



Q 



dQ 



dS,,, 



duo 



-Q, 



dGi 



dl 



—uo 
2 



(16) 



because Q-ball solutions give the extrema of Eq, S^, and Gi as regards Q, and /, respec- 
tively. These variables match the corresponding "thermodynamic" ones: Eq, uj, Q, Si^,a.ndGi 
correspond to the internal energy, chemical potential, particle number, and "thermodynamic" 
potentials [35]. After computing or Gj, one can calculate Q or ju"^ using the second or third 
relation in Eq. (fTBj) . and can compute Eq using Eq. (|TU|) or Eq. (fT^ . i.e. 

dS,.} 



Q 



duj 



Eq =ujQ + 



(17) 



or similarly Gj 



^ ^ Eq = Gi+ (icu2) /, S^ = Gi- (icu^) /. We shall make use of 



this powerful technique later. 



2.3 Q-ball equation and existence condition 

Let us consider the action S = — J dtSui in Eq. ([21) with our ansatz Eq. ([T3l) and the following 
boundary condition on a (Z^ — l)-sphere which represents spatial infinity 

a'\ = on the {D - l)-sphere. (18) 



Varying with respect to a we obtain the Q-ball equation: 

d'^a D — Ida dU,,, 



dr^ r dr da 



0, (19) 



d ( 1 fdaV , \ D-l fdaV , , 

There is a well known mechanical analogy for describing the Q-ball solution of Eq. ([T9ll [1] , and 
that comes from viewing Eq. ffTUl) in terms of the Newtonian dynamics of an unit-mass particle 
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with position a, moving in potential — f^j with a friction ^^f^, where r is interpreted as a time 
co-ordinate. Moreover pg = tua'^ can be considered as the angular momentum B Note that 
the friction term is proportional to ^^f^, hence becomes significant for high D and/or small r. 
According to Eq. ( l20l) . the "total energy", | (^)" — Ui^, is conserved for D = 1 and/or r — > oo, 
implying that in that limit the Q-balls have no radial pressure. Of course these are really 
field theory objects, and consequently, more restrictions apply: (i) no symmetry breaking, in 
other words a{r large) = 0; U"{a = 0) = > with an effective mass m, (ii) regularity 
condition: a'{r = 0) = 0, (iii) reflection symmetry under a —a. Note that Eq. fll9p coupled 
with the boundary condition Eq. (fTSl) implies a{r) is a monotonically decreasing function, i.e. 
a' < 0. In fact, according to Eqs. (fTSl [T9i) and the above conditions, our mechanical analogy 
implies that a static particle with a unit mass should be released somewhere on its potential, 
eventually reaching the origin at large (but finite) time and stopping there due to the presence 
of a position- and D- dependent friction. These requirements constrain the allowed forms of 
the U(l) potentials: for example if the local maximum of the effective potential — is less 
than 0, the "particle" can not reach the origin, a process known as undershooting. To avoid 
undershooting we require 

max{-U^) > ^ min < (21) 

If —Uuj is convex at a = 0, the "particle" cannot stop at the origin, a situation termed 
overshooting such that 



da^ 



<0 ^uj"^ < 



(22) 



(T=0 

Combining Eqs. (pTl [22]) . we find the condition on u for the existence of a single Q-ball: 



u!^ < \uj\ < (23) 



where we have defined the lower and upper bounds of u as ci^zp, i.e. uP^ = min [^) > and 
a;^ = ^^|o-=o = rn'^- The case, u;_ = 0, corresponds to degenerate vacua potentials (DVPs), 
while the other case, ^ 0, does not have degenerate vacua (NDVPs). The existence condition 
in Eq. fl^Hl) shows that U(l) potentials must have a non-linear interaction and is weakly 
attractive [3]. For convenience we define the maximum of the effective potential to be at 0"+ 
(i.e. ^^\a=a+ = 0), thus uj'i = and U^_{o-^) = 0. Moreover cr_ satisfies U^{cr_) = for 
o"_ 7^ 0. Notice (t_ ~ o"+ when u c:^ uj_. In Fig. [H we indicate the above introduced parameters, 
a±, uo_ using typical original and effective potentials for both DVP (left) and NDVP (right). 
To proceed with analytical arguments, we consider the two limiting values of u; or ctq = cr(0) 
which describe 

• thin wall Q-balls when ~ ct;_ or (jn ^ (7+ , , 
^ + (24) 

• thick wall Q-balls when uo ~ c<j+ or (Tq — cr_ . 

^/ is realised as an inertia moment in this mechanical analogy [U [5] • 
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We will not be considering Q-ball solutions that exist in a false vacua where u;^ < [9] or in flat 
potentials. When it comes to obtaining Q-ball profiles numerically, we will adopt a standard 
shooting method which fine-tunes the "initial positions" ctq subject to cr„(a;) < ctq < cr+('^); in 
order to avoid undershooting and overshooting. 
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B 



vacua potential (NDVP) with A = 1, B = ^ on the right. The DVP has degenerate vacua in 
line) where we set uj— = 0. The NDVP does not have degenerate vacua, but with lu = lu- = 
effective potential —Uui does have degenerate vacua. The two lines in the lower limit lu = ui- 
defined the maximum of the effective potential to be at (t+ and C/^j (p— ) = for a— ^ 0. The 
thick wall limit = 0;+ . With some values of (green dotted lines) satisfying the existence 
show the values of crq: clearly. 



a;+ = 1 and the effective potentials 
= I on the left and non-degenerate 
the original potential —IJ (red-solid 
= 0.5 (sky-blue dot-dashed line) the 
show that tT_ — ♦ cr+ where we have 
purple lines show a- — > with the 
condition Eq. H19|l. both potentials 



2.4 Three kinds of stability 

Absolute stability When the volume Vd approaches infinity [3] and/or lo is outside the 
limits of Eq. ( l23l) . then plane wave solutions can exist around the vacua of ?7(|0|). The equa- 
tion of motion for becomes a free Klein-Gordon equation whose solution can be written as 
= A^e*'-'''^"'^''*^ where uj^ = \/vn?- + and the normalisation factor = \J^^^ has been 
calculated from Q. Then the energy of the plane wave solution is proportional to ujk and Q 
linearly: Ef^ee = ^kQ -f'/ree — ''^Q where we have taken the infra-red limit, to obtain the 
second relation. The energy Ef^ee can be interpreted as the energy of a number Q of free 
particle quanta with the rest masses m. Furthermore, one might expect that the thick wall 
Q-ball energy with ^ u;+ = m approaches Efree because the Q-ball profiles approach zero 
exponentially at infinity [3]: 

Eq{uj = ^ Efree ^ mQ . (25) 
Hence the absolute stability condition for a Q-ball becomes 

Eq{u) < Efree. (26) 

We would expect Eq. fl^Uj) to be the strongest condition which a Q-ball solution has to satisfy. 
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Classical stability The classical stability [2], [3] can be defined in terms of tlie mass-squared 
of the fluctuations around a Q-ball solution. For zero mass fluctuations this corresponds to a 
zero mode, i.e. translation and phase transformation of the Q-ball solution. Using collective 
coordinates and Eq. (fT9l) which extremises S^^, such a mode can be removed. Since a detailed 
analysis can be found in the literature [21 [3], we simply state the final result which implies the 
classical stability condition is 

where we have used Eq. fll6l) in the second relation of Eq. (1271) . Since lo and Q have the same 
sign, the sign of ^ signals whether the solution is classically stable. The first relation of 
Eq. (1271) indicates the presence of an extreme charge in the parameter space of a;, (we will later 
see that the extreme charge at some critical value lo = turns out to be the minimum allowed 
and will be denoted by Qmin)- Let us remark on the characteristic slope of Eq/Q as a function 
of u: 

d ( Eq \ dQ 



duj\Q ) Q^duj-^ 

where we have used Eq. ( ITOj) and Eq. (fT6l) . Since, as we will see, > for all possible values 
of u within Eq. ( 123|) (in a specific potential [70]), the classically stable Q-balls should satisfy 
^ (^) — ^' "^^^ conditions from both Eq. ( 1271) and Eq. ( l28l) must be same. 

Stability against fission Suppose the total energy of two Q-balls is less than the energy of 
a single Q-ball carrying the same total charge. The single Q-ball naturally decays into two or 
more with some release of energy. As shown in [S], the stability condition against fission for a 
Q-ball is given by 

d'^En duo , , 

where we have used Eq. (1161) in the second relation of Eq. ( 129|) . Note that this is the same 
condition as we found above in Eqs. (ITTt [2HI) . so the condition for classical stability is identical 
to that of stability against fission. 

Trying to summarise the stability we can categorise three types of a Q-ball: i.e. absolutely 
stable, meta-stable, or unstable Q-balls. Absolutely stable Q-balls are stable quantum mechan- 
ically as well as classically; meta-stable Q-balls decay into free particle quanta, but are stable 
under small fluctuations; whereas unstable Q-balls sometimes called Q-clouds [70] decay into 
lower energy Q-balls or free particle quanta. 



2.5 Virial theorem 

Derrick's theorem restricts the existence of static non-trivial scalar field solutions in terms of 
a number of spatial dimensions. For example in a real scalar field theory, non-trivial solutions 
exist only in one-dimension, e.g. Klein-Gordon kink. Q-balls (or any non-topological solitons). 
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however, avoid this constraint because they are time- dependent (stationary) solutions [TTl [T2]. 
We can easily show this and in doing so obtain useful information about the scaling properties 
of the Q-balls as a function of dimensionality as well as the ratio between the surface and 
potential energies. Following [12], we begin by scaling the Q-ball ansatz Eq. (fT3l) using a one- 
parameter family r — > ar, whilst keeping Q fixed. Defining a surface energy S = Jy^ |cr'^) ^ 
potential energy U = jy U, and recalling that the charge satisfies Q = lu, we see that the 
energy of the Q-ball, Eq. ffTUj) . becomes 

Eq = S + U + ^. (30) 

dE' 

Now, under the scaling r — > ar, then Eq Eq where = because the Q-ball solutions 

a=l 

are the extrema (minima) of Eq. Evaluating this, we obtain the virial relation relating U and 
S 

DU = -{D-2)S + D^>0^Q'> g^n (31) 

where we have used our earlier notation, U > 0, for any values of a and defined the minimum 
charge (corresponding to W = 0) as Qmin = ^^^^-^5 > -v^ -D > 2. Since Q is taken as 
real and positive, no conditions appear for D = 1, 2. Notice this does not mean that one- 
and two- dimensional Q-balls do not exist, as can be seen from Eqs. (fT^ I^Hj) . The case of 
Q = recovers Derrick's theorem, showing no time- independent solutions for D > 2 [12]. 
Using S = ^{D-2 + D^y from Eq. (Ell), the characteristic slope Eq/luQ becomes 

l+(D-2 + D^y\ (32) 



cuQ \ S 

Let us consider three cases: (i)W ^ S, (ii)W — S, and (iii)W <^ S. They lead to predictions 
for the ratio of the Q-ball energy Eq to the energy contribution from the charge uQ: 



Eq 
uQ 



1 for 

2D~1 
2{D-1) 

§5i for (iii 



for (ii), (33) 



All of the Q-balls in the range of uj are classically stable because the terms, Eq/Q, monotonically 
increase as a function of cj, see Eqs. f l28l [H^ . The first case (i) corresponds to the extreme thin 
and thick wall limits a; ~ a;=p as will see. In the second case (ii), the potential energy is of the 
same order as the surface energy which means S and lA have equally virialised. This case will 
turn out to be that of the thin wall limit for DVPs when the surface effects are included. 

Suppose SjlA = const, over the large range of u within the existence condition Eq. (12^ 
except cij ~ where Eq/uj^Q ~ 1. We can find an approximate threshold value Ua for a 
Q-ball to be absolutely stable using Eqs. ( l25l ESj): 



1 for (i 

^ \ 2{D-1) 



m 



§5f for (iii). 



for (ii), (34) 
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Roughly speaking, Q-balls are classically and absolutely stable ifuKUa because of Eqs. ( l26l 1281) 
and Eq. (1331) . These approximations can and will be justified by our numerical results. We will 
find that the virial relation is a powerful tool enabling us to find appropriate values of Ua as 
opposed to the rather complicated computations we will have to adopt in the next section. We 
should point out a caveat in this argument, the assumption we are making here, that most of 
the Q-balls have an identical energy ratio S /U over a range of o;, does of course rely on the 
specific form of the potential. We have to remind the readers that the virial relation Eq. fl3ip 
gives only the relation between S and U, and states the presence of the minimum charge Qmin 
if the system allows the time-dependent (Q-ball) solutions Eq. ( fT3i) to exist. 



3 Thin and thick wall approximations 

In this section we will obtain approximate solutions for Q-balls in D-dimensions based on 
the well known thin and thick wall approximations for the radial profiles a{r) of the fields. 
Moreover we will show how we can then use these results to verify the solutions we obtained in 
the previous section for Eq/ujQ in Eq. flHHj) . Further we will then be able to test the solutions 
against detailed numerical solutions in section HI We start with two simple ansatze for the radial 
profiles, a step-like function for the thin wall case u; ~ tu- and a Gaussian function for the thick 
wall case uj ~ uj+. In both cases we will evaluate S"^, Q, Eq, as well as the conditions for 
classical and absolute stability before modifying the ansatze. Following that we will repeat the 
same calculations using our more physically motivated ansatze via the Legendre transformation 
technique described in Eq. ( IT71) . Let us comment briefiy on the form of the potential. We will 
see that in the thin wall limit ctq — cr+ with our modified ansatz, although in principle we do 
not have to restrict ourselves to particular potentials, we will not be able to investigate cases 
where the effective potential is extremely flat, hence we will have to limit our investigation to 
situations where this is not the case. In the thick wall limit uj ~ uj^, we have to restrict our 
analysis to the case of polynomial potentials of the form: 

U{a) = \m^a'^ - Aa'' + ^ B^a'' (35) 

p>n 

where n > 3, with the nonlinear couplings A > and -Bp > to ensure the existence of Q-ball 
solutions. We expect the thin wall approximation to be valid for general Q-ball potentials in 
which the Q-ball contains a lot of charge, with ~ cu^ > 0. In this limit we can define a 
positive infinitesimal parameter, e^^ = — f/(^((T+) ~ — ^-)^+ ^ O5 ^ind the effective mass 
around a "false" vacuum is given by, jj? = ^^lo-^. The other extreme case corresponds to the 
thick wall limit which is valid for Q-balls containing a small amount of charge, and it satisfies 
cij^ ^ cj^ = 711? . For later convenience, in this limit, we define a positive infinitesimal parameter, 
w?, = w? — Uj"^ > 
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3.1 Thin wall approximation for D > 2 

Step-like ansatz u; ~ At a first step, we review the standard results in the thin-wall 
approximation originally obtained by Coleman [1]. Adopting a step-like ansatz for the profile 
we write 

where Rq and ao will be defined in terms of the underlying parameters, by minimising the 
Q-ball energy. We can easily calculate S^j, Q, and Eq: 

= (Uo - lio'^al] Vo, Q = Lua^Vn, Eq = + UoVd (37) 

\ 2 J 2a^VD 

where Uq = U{o-q) and Vd = Vd^v = Rq). Note that Eq. (!37j) satisfies the Legendre transfor- 
mation results Eq. ffTTj) as we would have hoped. Since the ansatz Eq. (IHUjl neglects the surface 
effects, we are working in the regime U ^ S in Eq. 0331) . Therefore we should be able to re- 
produce the result, Eq ~ uQ with this solution. To see this, we note that the two terms in Eq 
are the contributions from the charge and potential energies.. These two contributions virialise 
since Eq is extremised with respect to Vd for a fixed charge Q, i.e. OEq/ OVdIq = 0, hence 

Vd = Q vV(2ap7o)- This then fixes Rq because we know, for a {D — l)-sphere, Vd = -^VLd~i 
where ^d-i = / (KId^i = r{D/2) ■ Substituting Vd = Q^/T/(2a^Uo) into Eq [the third equation 
of Eq. fl37|) ] and minimising Eq with respect to Cq, we obtain 

2 „2 




Eq = Q ■ min W = QuJ^ = uj_a^VD (3 



. Thus we recover Eq. ( 1331) in the 



where we have used Eq. ( I23l) in which cu^ = min y^j 

limit U ^ S. Finally, we remind the reader that we have obtained the minimised energy Eq 
with respect to Vd{Rq) and ctq in the extreme limit a; = a;_ where we find 



(To = a+. (39) 



Eq. (1391) implies that the "particle" spends a lot of "time" around cr+ because the effective 
potential — f/o; around cr+ is "fiat". Note that Q and Eq are proportional to the volume Vd in 
Eqs. f lH7|) and flHSl) just as they are for ordinary matter, in this case Coleman called it Q- matter 

pa- 

The modified ansatz ctq — (j^ Having seen the effect of an infinitely thin wall, it is natural to 
ask what happens if we allow for a more realistic case where the wall has a thickness associated 
with it? Modifying the previous step-like ansatz to include this possibility [301 H] will allow 
us to include surface effects [H [9], [10] and is applicable for a wider range of uj than in the 
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step-like case. Using the results, we will examine the two different types of potentials, DVPs 
and NDVPs, which lead to the different cases of Eq. flHHj) . 
Following [30], the modified ansatz is written as 

0-+ - s(r) for r < Rq, 
air) = { a(r) for Rq < r < Rq + 6, (40) 

ioT Rq + 6 < r, 

where as before the core size Rq, the core thickness 6, the core profile s(r), and the shell profile 
d"(r) will be obtained in terms of the underlying parameters by extremising S^^ in terms of 
a degree of freedom Rq. Continuity of the solution demands that we smoothly continue the 
profile at r = Rq, namely o-+ - s{Rq) = a-{RQ) and s'{Rq) = ^'{Rq). 

We expand U^^ to leading order around cr+, to give U^^^a) ~ — + |yU^s^ where s(r) = 
cr-|_ — o"(r). In terms of our mechanical analogy, the "particle" will stay around 0"+ for a long 
"time". Once it begins to roll off the top of the potential hill, the damping due to friction 
(oc [D — l)/r) becomes negligible and the "particle" quickly reaches the origin. Therefore, we 
can naturally assume 

Rq > S (41) 

where 5 is the core thickness. We know that cr'(O) = — s'(0) — O5 ^'(^q) O5 s'{r) > 0. 
Using Eq. ffTUl) . the core profile s(r) for r < Rq satisfies the Laplace equation: 

s" + ^^^s' -fi^s = (42) 
r 

whose solution is 

s{r) = r(i-f ) (Ci/D_,(/ir) + C2A'D_,(/ir)) (43) 

where / and K are, respectively, growing and decaying Bessel functions, Ci and C2 are con- 
stants. Since s(0) is finite and s'(r) > 0, it implies that C2 := 0. Since Iu{z) ~ /2V{i' + 1) 
for small z = fir and u —1, —2, —3 . . . , thus s(0) is finite: 

which gives a relation between Ci and ctq. Also the analytic solution is regular at r = 0: 
s'(0) ~ 0. For large r ~ Rq, Eq. (gSD leads to 

5' ^-2 ,,,, 
/i . /i. (45) 



s 



where we are assuming 



/i > 1/Rq (46) 



12 



and have used the approximation Iu{z) ~ for large z = nr. As already mentioned we 
note that this result is not strictly valid for extremely flat potentials, i.e. /i ~ 1/Rq, because 
the expansion is only valid for z = /ir ^ 1. We will therefore only be applying it to the cases 
where the effective potential is not very fiat. 

Turning our attention to the shell regime Rq < r < Rq + S. Considering the "friction" 
term in Eq. (fT9i) . we see that it becomes less important for large r compared to the first and 
third terms in Eq. ( 142|) . because 



D-1 



Rc 



s'{Rq) 



D 



flRq 



■fi^siRq) 



« fi^siRg) ^ s"{Rq) 



dU.., 



ds 



(47) 



--Rq 



where we have made use of Eqs. ( l45| 1461) . Imposing continuity conditions, namely 
ai^Rq), —s'{Rq) = a'^Rq), Eq. ( fT9i) without the "friction" term becomes 



d'^a dU„: 



dr"^ da 



0, 



■s{Rq) 



(4^ 



where a{r) is defined as being the solution to Eq. fHHl) . With the condition ^{Rq) = a^ — s{RQ) 
and Eq. dS]), we find aiRg) = (t+ - y^^^^^^^ffD^iw^ • For the thin wall limit (Tq ~ a+, 
we obtain cr(i?Q) ~ a+ — A where A is an infinitesimal parameter which satisfies cr+ ^ A. 
Therefore 

a{RQ) » s{Rq). (49) 

Although Eq. does not hold exactly, the "total energy", ^ (^)^ - ~ with Eq. (fT8|), 
is effectively conserved with the radial pressure pr vanishing outside the Q-ball core. This 
fact implies that the surface and effective potential energies virialise with equal contributions, 
Ssheii — Usheii — shell where we have introduced shell and core regimes defined by Xcore = 
VLd-i jl^^ drr^~^ ■ ■ ■ and Xgheii = ^D-i Irq^^ drr^~^ ■ ■ ■ for some quantity X. Using a' < 
and the condition d'{RQ + 6) = 0, the thickness of the Q-ball can be written as 6{uj) = 



^{Rq) da 



■■. Since 6 is real and positive, we have to impose 



(t{Rq) < o-_. 



(50) 



recalling ?7(^(cr_) = for (t_ 7^ 0. 

With the use of Eq. (|T7|) . we turn our attention to extremising the Euclidean action 
in Eq. fill I) for the degree of freedom Rq. Using the obtained value Rq, we will differentiate 
5*^ with respect to uj to obtain Q as in Eq. (fT6l) which leads us to the Q-ball energy Eq as 
in Eq. ( |TOl) and the characteristic slope Eq/ujQ. For convenience we split 5*^^ into the core 
part S^°^''^ for r < Rq and the shell part S^^'^'"'' for Rq < r < Rq + 6 using Eq. ( HOl) . Using 

Vd = ^VLd-i > dVn = R^-^flD-i > d^Vo = R^-^VLd-2 and Eqs. dH SS]), we find, 

= -Vn ■ 6. + dVn ■ (i/is^(i?Q)) - d'Vo ■ (^^^^^/^^'(^q)) (51) 
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where the first term, e^^, in Eq. ( l5Ti) comes from the effective potential energy, while the second 
and third terms arise from the surface energy. Since is an infinitesimal parameter in the 
other thin wall limit ^ cu^, it gives 

The effective potential energy balances the surface energy in the shell (see Eq. fl48p ). therefore 
by introducing the definition T = Jq^^'^^ dayj2U^, we see 

S'f'^" = fiz)-i / dar''-^^2Uja)<nD-i{RQ + 5)''-^T 
Jo 



(53) 



dVD-T + d'Vo ■ ( ^^{D -1)6.t)+ 0{Rg-\ ■ T, (54) 



^(Z.-l)..T)+0(i.-| 

where we have used the fact that the integrand has a peak at r = Rq + 5 in the second relation 
of Eq. (l53l) [72] and Taylor-expanded {Rq + 6)^~^ in going from Eq. ( l53i) to Eq. ( l54l) because 
of our approximation Eq. fl4ip . Combining both expressions Eqs. 051^ [3^ . we obtain 

5'. = S'r^ + S':'^'^" (55) 
^ -e^-VD + T-dVo + h-d^Vo (56) 

where r = T + ^^s'^{Rq). Note that whilst in r, T contains the equally virialised surface and 
effective potential energies from the shell, the second term ^^s'^{Rq) contains a surface energy 

term from the core. Moreover we have defined h = 



n 



D-l 



{D - 1)6 ■ T - ^-^I^s\Rq) 



o \^ ^ which 

is negligible compared to r because of the assumptions Eqs. (HTl H6l) . Therefore, we will take 
into account only the first two terms in S^j, Eq. ( 1561) . It is also important to realise that 



/ da^/2U^+ / da^2U^_ / da^2U^_ = const. (57) 

Jo Ja(Ro) Jo 



which is independent of u and D in the limit of — > a;., where we have used the other thin 
wall limit u; ~ Our modified ansatz is not only valid in the extreme limit u = but also 
in the limit ~ as long as r depends on u "weakly". Note that the condition of Eq. (l50l) 
also ensures that r is positive and real. In addition, the second term in Eq. (1571) is negligible 
compared to the first term, i.e. 

Sshell — Ushell — -^^Qshell ^ Score (58) 

because of cr+ ~ o"(-Rq), see Eq. ( 149|) . 

We can make progress by using the Legendre transformation of Eq. (JTTI) . which implies 
that we need to find the extrema of S'j^ with fixed u, i.e. = 0. This is equivalent to the 
virialsation between e^^ and r. Then one can compute the core radius, 

Rq = (D-1)-. (59) 
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Note that this imphes that one- dimensional thin wall Q-balls do not exist due to the positivity 
of Rq. By using Eqs. fl3Ul EH]) and Eq. ffT7|) . we can compute the desired quantities to compare 
with the results we obtained using the step-like ansatz, in particular Eqs. ( 1371 [38|) . and we can 
confirm that the classical stability condition Eq. ( 1271) is satisfied: 

~ LdVD = ^^VD>0, QH-^a^V^ (60) 

Eq ~ uj^alVD + ^dVo (61) 



2D -I 



2(D-1) 2(D-l)a;2 



(62) 



--^ ~ 1 ^ ~ i < 0. (63) 

Q dto €^ €^ 

We can see the virialisation between r and e^^ for the second and third terms in Eq. ( l60l) . As in 
Eq. ( l38l) . the first term of Eq, in Eq. ( 16T1) . is a combination of an energy from the charge and 
potential energy from the core throughout the volume, while the new second term called the 
surface tension, represents the equally virialised surface and effective potential energies from the 
shell as in Eq. (|58l) . In the limit u; ~ e^^ becomes asymptotically zero which implies Eq. ( !52l) 

as we saw. We have also seen S shell > Score- Using U = Ucore+UshelU <S = Score + <S shell ~ S shell, 

and Eqs. (!52ll58l). we obtain 

U r^S + uj_Q (64) 

which we will use shortly. Since the characteristic function, Eq/Q, increases monotonically as 
a function of uo and S*^ > 0, i.e. ^ {~Q^ > or Eq. (lUHj) . the classical stability condition 
Eqs. (|271 [28]) is satisfied without specifying any detailed potential forms. However, the physical 
properties of the finite thickness thin wall Q-balls do depend on the vacuum structures of the 
underlying potential. To demonstrate this we consider two cases of non-degenerate vacuum 
potentials (NDVPs) with cu- 7^ and degenerate vacuum potentials (DVPs) with = (see 
red-solid lines in Fig. [T]). Suppose that the thin wall Q-balls have identical features over a 
large range of cu, we can find the approximate threshold frequency uOa using Eqs. (^31 [H^ as we 
assumed when we obtained Eq. 



NDVPs This type of potential reproduces the results we obtained in Eq. (I38p corresponding 
to the regime U ^ S which corresponds to the existence of Q-matter in that the charge and 
energy are proportional to the volume Vd due to the negligible surface tension in Eq. (16T1) . 
Hence, the modified ansatz Eq. ( HOl) can be simplified into the original step-like ansatz Eq. ( 1371) 
with negligible surface effects in the extreme limit uj = uj^. To see that, we need to recall the 
definition of uj_ in Eq. ( 1211) . We can realise that /i is the same order as u;_ except the case of 
uj_ = 0. Using /i ~ we can show that ^uQ ^ Score ~ ^fJ'S'^{RQ)dVD where we have used 
Eqs. (l46l l49l). Using Eqs. (152 | l58l) and ^uQ ^ Score which we just showed, we can obtain the 
desired result U ^ S. Similarly Eq. ( l62l) in the limit u ^ lu^ simplifies to give ^ ~ 1 which 
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is the result of Eq. ( l33l) with U ^ S. Using Eq. (l62l) and Eqs. (125| 1331) . we can also find the 
critical value Ua for absolute stability 



D-l L {2D - 1) 

Finally, thin wall Q-balls in NDVPs are classically stable without the need for the detailed 

potential forms, however the absolute stability condition for a; ~ a;_ depends on the spatial 
dimensions D and on the mass m. 

DVPs For the case of the presence of degenerate minima where U- = 0, since e^^ = ^u'^a^, 
we immediately see from Eq. (I6T]) that 

^" '''-\ (66) 



uQ 2{D - 1) 

which reproduces Eq. flH^ for the case of 5 ~ W. As in NDVPs, we know Eq. flMl) in the 
limit u; ^ cij-, but the second term uj_Q becomes zero in the present potentials. It follows 
that Ucore — and Usheii — Ssheii ^ Scores heucc S hi. In other words, most of the Q-ball 
energy is stored within the shell. In addition the charge Q and energy Eq are not scaled 
by the volume, which implies that the modified ansatz does not recover the simple ansatz as 

opposed to NDVPs. In particular we find that Q = Rq~^^'^^ oc i?[f Eq = 

^^D-iT Rq-^ oc R^-\ a nice check of our general results follows by writing Eq in terms of 
the charge Q by eliminating Rq between the two expressions. This gives Eq oc Q'^(^~^y('^^~^) 
which reproduces the three dimensional results obtained in [30] . 

Finally, let us recap, the key approximations and conditions we have made in this modified 
ansatz. They are Eqs. (HT| HUl [30]) . and Eq. (IFTl) . We believe that the estimates we have 
arrived at for the thin wall Q'balls are valid as long as: the core size is much larger than the 
shell thickness; the effective potential is not too flat around cr+; the core thickness 6 and surface 
tension t/D are positive and real; r is insensitive to both u and D. With the extreme limit 
— > the Q-balls in DVPs recover the simple step-like ansatz, while the ones in NDVPs do 
not. One-dimensional Q-balls do not support thin wall approximation due to the absence of 
the friction term in Eq. (fT9|) . 

3.2 Thick wall approximation for u'^ ~ cj^, i.e. 

Gaussian ansatz As we have started with the simple step-like ansatz in the thin wall ap- 
proximation, a Gaussian function is a simple approximate proflle to describe the thick wall 
Q-balls [I2]. Using a Gaussian ansatz 

(T(r) = (Toexp ^--^^ , (67) 
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we will extremise S*^ with respect to (Tq and R with fixed cu, instead of minimising Eq with 
fixed Q. Notice that the slope —a' jo becomes 2rjR? which is linearly proportional to r and 
the solution is regular at r = 0: a"'(0) = 0. Neglecting higher order terms Bp in Eq. ( 135|) with 
Eq. ( 1671) . one can obtain straightforwardly 



Q 



TX 



D/2 



D / 2 

-m,^, H - A(Jn~^ 



n 



Er 



D/2 



Q 

UJ 



(68) 
(69) 



Eq. (fTTj) can be easily checked in Eqs. (l68l [69]) . The first and last terms in Eq. (169!) are the 
potential energy terms. The second energy term comes from the charge energy, and the surface 
energy term appears in the third term. By finding the extrema of S^^ with respect to ctq where 



= 0, it defines the underlying parameter ao as 



ml + 



2D 



1 

nA 



2 



l/n-2 



\2A 



l/n-2 



(70) 



where we have neglected the surface term and used the approximation D/2 ~ 0{1) in the 
second relation of Eq. ( 170|) . Then we can check that the Gaussian ansatz naturally satisfies the 
other thick wall limit ctq — o"- and that the higher order terms Bp in Eq. ( 135|) are negligible. 
Using the first relation of Eq. (170 p . one needs to extremise S^^ with respect to another degree 

which determines R: 



of freedom R with 



R 



l2{2-D) 



> 0. 



(71) 



The reality condition on R implies that the Gaussian ansatz is valid only for D = 1. The width 
of the gaussian function R in Eq. flTTl) becomes very large in the thick wall limit — > 0, thus 
we can justify that the surface term in Eqs. flU^ [70]) are negligible. Therefore we are looking 
at the regime U ^ S which should lead us to Eq ~ uoQ as in the first case of Eq. ( l33l) . To do 
this for D = 1 we substitute Eq. ( 1701) into Q, Eq, S^^: 



Q 

Eq 



1 1 

2^n 



1 

n ^ 
2m2 



2mlQ 



UJ 



> 0, 



1 



(72) 
(73) 



where we have considered the thick wall limit u; ~ m in the second relation of Eq. flTHl) . We can 
check Eq. ( !33l) and the analytic continuation Eq. ( l25l) . In the same limit, the Euclidian action 
becomes an infinitesimally small positive value: S^j ^ 0~^- 



17 



Using the second relation ctq in Eq. ( ITOl) and Eq. ( 17T]) . one can find 

Qduo ml\n-2 ) ml\n-2 )- ^ ' 

where we have used the fact that m^; is a positive infinitesimal parameter in the hmit, u; ~ Ci;+ 
going from the first relation to the second one. Eq. ([71]) shows that the classical stability 
condition clearly depends on the non-linear power n in the potential Eq. flH3|) : n < 6. This 
is contradictory because Eq. flTHl) gives ^ ("^) ^ ~^ ~'~ n '^hich implies n < 4 for the other 
classical stability condition using Eq. fl28l) . We will shortly see that this contradiction between 
Eq. (1271) and Eq. fl^Hj) is an artefact of the Gaussian ansatz. Moreover, our conclusion should 
state that the Gaussian approximation is valid only for D = 1. These awkward consequences 
are improved with the following physically motivated ansatz. 



The modified ansatz Having considered the case of the simple Gaussian ansatz following 
the spirit of [I2], we found some problems for the classical stability. To fix these, we need a 
more realistic ansatz [TTl [30], [SS] [39] [52]. To do this we drop an explicit detailed profile to 
describe thick wall Q-balls and rescale the field profile so as to work in dimensionless units 
whilst extracting out the explicit dependence on u from the integral in Si_j. As in the thin wall 
approximation with the modified ansatz, we will again make use of the technique Eq. (fT7j) to 
obtain other physical quantities from S^. 

We begin by defining a = aa and r = br with a and b which will depend on u. Substituting 
them into Eq. ffTT]) with the potential Eq. flH^ we obtain: 

= |rfl],,„i|rfrT^-^6^|i(^)V'2 + ia2m^a^-Aa"a"+ (75) 



p>n 



i^yno-i I dfr^~^]^ + ^2 - a" + Y,Bpb''aF-''~aA , (76) 



m 



4/(„.2)-D+2^-2/(n~2)^^_^^^ (77) 



with the rescaled action Sn = f drr^ ^ {\^'^ ^) ^'^^^ ^ ~ ^ i*^"' ^^'^ have neglected 
the higher order terms involving Bp. In going from Eq. (l75l) to Eq. ( 176]) we have set the 
coefficients of the first three terms in the brackets to be unity in order to explicitly remove the 
u) dependence from the integral in S^^. In other words we have set | (|) = \o?m'^ = Aa"'. This 

implies, a = ( ) = and b = . Crucially Sn is independent of on, and is positive 

definite [TTl [521 [30] . Adopting the powerful approach developed in Eq. ([T7]) . given S^^ we can 
differentiate it to obtain Q and then use the Legendre transformation to obtain Eq. This is 
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straightforward and yields 
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toQ. 



(78) 
(79) 

(80) 



The first term in Eq. (IHUj) is the energy contributed by the charge, while the second term is 
dominated by the effective potential energy, hence U ^ S. Therefore, we can also recover the 
result Eq ~ ujQ in the thick wall limit ~ as we would expect from Eq. (l33ll when U ^ S. 
Since Q and Eq should be positive definite, it places the constraint [52] 



D < 



n-2 



+ 2. 



With the condition Eq. (IHTj) . it is easy to see that Si^ in the thick wall limit, u; ~ cj^ 

where — 0"^. There is another constraint emerging from the need for the solution to be 
classically stable: 



UJ dQ ^ 
Q duj 



D 



n 



ijj 



^ D < 



n-2 
4 



D] <0. 



n 



12) 
i3) 



which coincides with Eq. (1741) in the case oi D = 1. Notice that the modified ansatz is valid 
not only for D = 1 but also D < + 2 in Eq. flHTl) . For D = 3 this result matches that 
of [Snj. The classical stability condition, Eq. flH^ . is consistent with the need for Q and Eq 
to be finite. Eq. ( !83l) is more restrictive than that given in Eq. ( IHTl) . Furthermore, we should 
check the relation Eq. ( l28i) for the characteristic function Eq/Q in terms of u. It follows that 
^ (^^) ~ 1 — 2 — D + 2) ^ > 0, which requires the same condition as Eq. flH^ . With 
this fact and Eq. flHOjl . it implies that thick wall Q-balls satisfy the conditions for both classical 
and absolute stability. Moreover it also reproduces the previous results in [73], for the case 
of -D = 2 and n = 4, p = 6 (6-th order potential). Unlike the Gaussian ansatz Eq. flU7|) . our 
modified ansatz now shows consistent results between Eq. (127|) and Eq. ( l28l) . 

Let us remark on the validity of our analysis following [52]. In this section we have used a 
modified ansatz which has involved a re-scaling of a and r in such a way as to leave us with 
a dimensionless action S'„. There are restrictions on our ability to do this as first pointed out 
in [52] for the case of D = 3. We can generalise this to our D dimensional case. Given that 
the Q-ball solutions extremise Sn, we may rescale r or cr introducing a one-parameter rescaling, 
r ^ ar OT a ^ Xa which will deform the original solution. Defining X{a) = Sn[o:r, a{ar)] and 
Y{X) = Sn[Xa{r)], we impose the condition that the action Sn is extremised when a = A = 1, 
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which imphes = = ^\x=i. It is possible to show that these conditions imply that 

consistent solutions require the same condition as Eq. flHT]) . The three dimensional case leads 
to the result, n < 6, as originally obtained in [30]. The particular choice of n = 4 which we 
will investigate shortly implies D < A for the validity of our thick wall approximation with the 
modified ansatz. Moreover, thick wall Q-balls become classically unstable for D > 3 a.s can be 
seen from Eq. ( !HH|) . 

What have we learnt from extending the ansatz beyond the Gaussian one? We have seen 
that they have lead to different results. For instance, the Gaussian ansatz essentially has a 
contradiction for the classical stability analysis even for D = 1, whereas the solutions based 
on the modified ansatz are valid for D, which satisfies Eq. flHTj) . and give consistent results, 
Eq. f lHHl) . for classical stability. 

As we will see in the next section, our numerical results in which we obtain the full Q-ball 
solution support the modified ansatze for both thin and thick wall cases. 



4 Numerical results 

In this section we obtain numerical solutions for Q-balls using the polynomial potential in 
Eq. (1351) including only one higher order term: 

U{a) = ^m^a^ - Aa" + Ba^, (84) 

where A > 0, -B > 0,p > n > 2. We shall confirm the results obtained analytically using the 
modified ansatze for both the thin and thick wall Q-balls. Recall that Uuj{o') = U{a) — ^u'^a^, 
with f/^((T_) = and cr+ marks the maximum of the effective potential — f/^ where cr_|_ 7^ 0. 
For a particular case p = 2{n — 1) , we find 

, \ l/(n-2) / , \ l/("-2) 

c^-M=| 1^ J .o^{-)=y I • (85) 

Also, for convenience, we set 



cu+ = m = l, cj„ = Wl- — >0^A< V25, (86) 
V 2x5 

where we recall the definitions of u;+ and u;_ are that u"^ = ^^\^=q = m? and f/tj_(c'"+) = 0. 
Setting cu. = in Eq. ( l86i) implies that U{a) in Eq. (l84l) has degenerate vacua at a = 0, ±cr+, 
whilst the original potential U does not have degenerate vacua with uj^ 7^ 0. In this section, 
we shall consider two examples of the potential f/, which can be seen as the red (solid) lines 
in Fig. [TJ The degenerate vacua potential (DVP) on the left has a;_ = (A = V2B) and the 
non-degenerate vacua potential (NDVP) on the right has uj- = 0.5 {A = \j3B/2). In order 
to determine actual values for A and B, we define cr+(a;+) = 1 and set = 4, p = 6 for both 
cases, hence A = |, i? = | in DVP and A = 1, i? = | in NDVP. Fig. [T]also includes plots of 
the effective potentials for various values of uj. 
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4.1 Numerical techniques and parameters 

To obtain the Q-ball profile we need to know the initial "position" ao = o"(r = 0). This is done 
using a shooting method, whereby we initially guess at a value of (Tq, then solve Eq. (fT9l) for the 
Q-ball profile, and depending on whether we overshoot or undershoot the required final value 
of (T, we modify our guess for ctq and try again. Throughout our simulations, we need to specify 
the following three small parameters, e, C,, r] which, respectively, determine our simulation size, 
rmax, the radius at which we can match the analytic and numerical solutions (Rana), and the 
core size Rq. The smoothly continued profile is computed up to r = i?„ 



''max ■ 



Shooting method Let us consider an effective potential — f4j which satisfies the Q-ball 
existence condition, Eq. ( 1231) . We have to initially guess ctq subject to it be being in the 



appropriate region a- < ao < a+. For example it might be = . There are then 

three possibilities, the particle could overshoot, undershoot, or shoot properly. The last case is 
unlikely unless we are really "lucky". If it overshoots then we would find o"(ro) < at some 
"time" To- If that were to happen we could update to = "g^^- g^g q^j. j^g^t guess. On the 
other hand if it undershoots, the "velocity" of the "particle" might be positive at some "time" 
r(7, cj'(r[/) > 0. If that were to happen we might update a% to o\< = "^^^g q^j, v^&^i guess. 
After repeating the same procedures say times, we obtain the finely-tuned initial "position" 
cTg — 0"^ as our true value. To be compatible with numerical errors, our numerical simulation 
should be stopped with an appropriate accuracy parameterised by e: 

e > a{ru = Vmax) > (87) 

where Vmax is the size of our simulations, and e measures the numerical accuracy where a small 
value of e corresponds to good numerical accuracy. Unfortunately the final profiles still have 
small numerical errors for large r. To compensate for these errors, the profiles should continue 
to the analytical ones smoothly at some point r = Rana- 



Matching analytic and numerical solutions at Rana For large r, the Q-ball Eq. (HM can 

be reduced to 

a" + ^^-^(y' -mla = Q. (88) 

The analytic solution becomes 



where E is a, constant which is determined later. Note that we have used the fact that the 
modified Bessel function of the second kind has the relation K^{r) ~ V^e"'' for large r and any 
real number /i. In order to smoothly continue to the analytic profile Eq. ( 189|) at the continuing 
point Ranai the following condition is required using the second relation of Eq. (!89|) : 



D-l aL 



, ^ I num 

+ + 



2r cr„ 



< e (90) 
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where a parameter ^ should be relatively small. Hence we can find the appropriate profile in 
the whole space 



for r < Rana, 

a (H ] ( Roiml)^^"-^^/'^ f-"^'^{r-Rana) fori? <r<R 



(91) 



where we have computed E using Eq. flHOl) and our simulations are carried out up to r = 



Core size and thickness of thin- wall Q-ball Using Eq. we can define the core size 
r = Rq and the numerical thickness 5num. by the slope —a' jo with the following condition 



D 



G 



cr-n 



(92) 
(93) 



Notice that the definition of 5num is different from the definition in Eq. (jUj) where 

(^(Rq) da 

Jo 



Numerical parameters We have run our code in two different regimes of uj for both DVP 
and NDVP because the profiles for large uj are needed to look into larger simulation size r^ax 
compared to the ones for small uj. Due to numerical complications, we do not conduct our 
simulations near the extreme thin wall limit, i.e. ~ lj^. However, by solving close to the thin 
wall limit, our numerical results for ctq — o-+ and -Rq ^ 5num allow us to recover the expected 
properties of thin wall Q-balls with the modified ansatz Eq. fjlUj) . Finally, our results presented 
here correspond to the particular sets of parameters summarised in Table: [H 



DVP 


UJ 


e 


^max 


Rmax 




V 


0.38-0.73 


4.0 xlO-2 


30 


200 


8.0 xlO-3 


1.0 xlO^i 


0.73-0.99999 


1.0x10-5 


40 


200 


8.0 xlO-3 


1.0 xlQ-^ 



NDVP 



UJ 


e 


'^max 


Rmax 




V 


0.60-0.85 


3.0 xl0~3 


30 


200 


8.0 xlO-3 


1.0 xlO-^ 


0.85-0.99999 


1.0x10-5 


50 


200 


8.0 xlO-3 


1.0 xlQ-^ 



Table 1: The numerical parameters in DVP (top) and in NDVP (bottom). 



4.2 Stationary properties in DVP and NDVP 

We devote a large part of this section to justifying the previously obtained analytical results in 
the thin and thick wall approximations by obtaining the appropriate numerical solutions. 
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Profiles with our numerical algorithm In the top two panels of Fig. [2] the two red lines 
(one dotted and one with circles) show the numerical slopes —o' jo for the case of i5 = 3 for 
two values of uj. These are then matched to the analytic profiles (green dotted lines) in order 
to achieve the full profile as given in Eq. ( l9Ti) . Recall that we expect in general for all values of 
(x;, the analytic fits to be accurate for large r, the numerical fits to be most accurate for small 
r and there to be an overlap region where they are both consistent with each other as seen 
in Fig. [21 We have also plotted in dot-dashed purple lines our analytic fits, Eq. (1^ . for the 
slopes of the thin wall cores from r = 0.5. We should remind the reader that this fit only really 
works for the case of small lo because we are dealing with thin wall Q-balls. Notice, it is clear 
from the purple lines that the core sizes cannot be determined by this technique for the case 
c<j = 0.9 ~ 

The bottom two panels show the full profiles satisfying Eq. ( jOTll for arbitrary D up to D = 5. 
We have been able to obtain the Q-ball profiles in the whole parameter space lo except for the 
extreme thin wall region ~ Both DVP and NDVP Q-balls have profiles with similar 
behaviours in that as the spatial dimension increases, so does their core size. 
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Figure 2; The top two panels show the numerical slopes ~a' ja for the case of D = 3 for two values of uj for both DVP (left) 
and NDVP (right). The red (one-dotted and one with circles) lines show the numerical slopes and the green dotted lines with two 
different widths the corresponding analytic solutions. The purple dot-dashed lines with two different widths show the analytic fits 
for the core profiles. The bottom two panels show the full Q-ball profile as described in Eq. H91| l for a number of values of a; and 
D. Note how the core size increases with D. 
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Criteria for the existence of a thin wall Q-ball with core size Rq The top and middle 
panels of Fig. E] show the numerical results for cro(ti;) and Snum/Rq against uj for a number of 
spatial dimensions D. For the case of il* > 3 it is clear from the top panels that the Q-balls are 
well described by the thin wall result Eq. (l85l) for most values of cu, with the range increasing 
as D increases. The case of Z) = 2 is less clear, it appears to asymptote onto the line. We 
believe there is a solution that exists for that case for small values of uj. An important point is 
that for the approximation to be valid we are working in the regime ^numlRq < 1 which can 
be seen to be true from the middle panels (again we believe the case of -D = 2 is heading below 
the line Snum/RQ = 1 for small uj. 

These results are consistent with our analytic solutions for finite thickness Q-balls given by 
Eq. fHUl) . subject to the criteria ctq — (7+ and Rq > Snum, even though u; ~ 

For D = 1 we see in the top panels that ao exactly matches cr_, (the orange dot-dashed 
lines). The bottom two panels in Fig. [3] show the core sizes Rq of thin- wall Q-balls which satisfy 
our criterion Eq. ( l92i) . Recall that Rq in Eq. (159|) is a function of uj assuming r depends on u 
weakly, thus we plot the numerical core sizes comparing them with our analytical approximation 
for DVP and NDVP, respectively 



2(-D l)TVium_ j^NDVP ^ 2(Z) 1)t, 

' ^ (u;2 — uj 



where the parameter Tnum is computed numerically (see Table: The presented numerical core 
sizes match excellently with the analytical fittings over a wide range of uj. Some numerical errors 
appear around uj ~ uj+ since we cannot determine the thick wall cores with this technique, see 
the top two panels in Fig. O Table: [2] shows analytical and numerical values of r using Eq. ( 1571) 
and the above fitting technique. We confirm that the values of r (a part of the surface tension 
t/D in Eq. flUTj) ) are nearly constant, depending slightly on D. Therefore, the assumptions we 
made for thin wall Q-balls are valid as long as ctq ~ cr+ and Rq > Snum- 

Configurations Fig. H] illustrates the configurations of charge density pg (top) and energy 
density pe (bottom), in both DVP (left) and NDVP (right). Each of the DVP energy densities 
around uj ^ uj- has a spike within the shells, while those spikes are not present in NDVP. The 
presence of spikes can contribute to the increase in surface energy S, which accounts for the 
different observed ratio for S /U in the two cases, where U is the potential energies. Otherwise 
DVP and NDVP models have similar profiles in Fig. [21 Moreover, we have numerically checked 
that Q-balls for D > 2 generally have positive radial pressures, whereas the ID radial pressures 
are always zero, i.e. ^a'"^ = U^j due to Eq. ( lT5l) . 



Virialisation and characteristic slope Eq/ujQ The top panels in Fig. El illustrate the 
ratios S /U and the four bottom ones show the characteristic slopes of Eq/ujQ against uj in 
both the thin-wall (middle-panels) and thick-wall (bottom-panels) limits. According to our 
analytic arguments Eq. (16^ . we expect S /lA 1 in the extreme limit ~ = in DVP. 
Similarly, we expect S /U ~ in the same extreme thin wall limit uj = uj_ = 0.5 for NDVP. The 
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Fig UrG 3i The initial "positions" (Tq (top), Snum/RQ (middle), and the core sizes Rq{uj) (bottom). The top panels show cr±, 
Eq. II85I I as black and orange dot-dashed lines respectively. The middle panels show the range of values of u> for a given value of D in 
which the core thickness is smaller than the core size, a crucial assumption we have to make. In the bottom panel, the analytical core 
sizes in Eo. l(94)l are plotted with the numerical ones for the following uj ranges: [0.38-0.40], [0.38-0.55], [0.38-0.60], [0.38-0.70] 
in DVP, and [0.60 - 0.62], [0.60 - 0.65], [0.60 - 0.75], [0.60 - 0.85] in NDVP and for D = 2, 3, 4, 5, respectively As can be seen, 
the fits are excellent. The range of uj values chosen have been based on the results shown in the top two panels and correspond to 
that range where the thin- wall Q-balls are solutions (except for D = 2). 

latter case corresponds to the existence of Q-matter with the simple step-like ansatz Eq. 
Although we are unable to probe these precise regimes, we believe the slopes of the curves 
indicate they are heading in the right direction. The thin wall slopes Eq / ujQ in the two middle 
panels lie nearby the analytical ones, Eqs. (I62l 1661) . as long as ctq — cr^ (see Fig. [3]) except for the 
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Table 2 ; The values of Tana and T„um in terms of D in DVP and NDVP, see Eqs. I I57I I59I I. 
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Fi gUre 4 1 The configurations for charge density pQ (top) and energy density pE (bottom) computed using Eq. I I14I I for both 
DVP (left) and NDVP (right). The presence of spikes of pE in DVPs contributes to their increased surface energies. 



2D cases because for D < 2 the profiles are not well fitted by thin wall predictions. Similarly, 
the thick wall slopes Eq/ujQ in the bottom two panels agree with our analytical predictions 
Eq. f l80l) using the modified ansatz rather than with Eq. flTSl) using the simple Gaussian ansatz. 
We have confirmed that the analytic thick wall slopes with Eq. flHUl) can not apply to higher 
dimensions D > 4, see Eq. ( IHTi) . Around the thick wall limit uj ~ uj+, the behaviours in both 
potentials are S <^IA (see top panel), which implies Eq ~ ujQ as predicted in Eqs. (173 1 179]) ; 
hence we can verify that the solutions are continued to the free particle solutions, see Eq. (!25|) . 
Our physically motivated modified ansatze in both the thin and thick wall limits therefore have 
clear advantages over the simple ansatze in Eqs. ( j36l 1671) . 

Q-ball stability Fig. [H] shows the classical and absolute stability lines for Q-balls. Table: [S] in- 
dicates the approximate analytical values of Ua derived by Eqs. (IMl 1651) . which can be compared 
to the numerically obtained critical values oj for the stabilities denoted by c^c, ojg, uJch, ^a, and 
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Fig U.r6 5 '. The ratio of S /lA where S and lA are surface and potential energies (top panels) , the characteristic slope Eq /ijjQ in 
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panels), with the analytic lines Eqs. H73I I80|l . 
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Uf in Table: HI Each of these are defined by 



duj I 





_ d (Eq\ 







0, Eq/Q\^^ 



^ch 



= respectively. The 3D analytical plots of ^ (^) in the thin and thick wall limits, 

Eqs. ([HSl IHHj) . can be seen to match the corresponding numerical data in the appropriate limits 
of uj. We have confirmed numerically that for both DVP and NDVP cases = i-^/ — "-^s — ^ch, 
see Table: HI This can be easily understood from Eqs. (1271 1281) and Eq. ( l29l) . 

Recall Eq. (l83l) . the polynomial potential Eq. (18^ with n = 4, leads to the classical stability 
condition D < 2 for the thick wall case. The top panels in Fig. E] demonstrate that thick wall 
Q-balls in D > 3 are classically unstable. In Table: HI one can check that the absolute stability 
condition is more severe than the classical one, then there are three types of Q-ball [2] as 
before: absolutely stable Q-balls for u < uja-, meta-stable Q-balls for oJa < oj < Uc, which are 
not quantum-mechanically stable but classically stable, and can decay into multiple Q-balls, or 
unstable Q-balls for uj^ < uj. 
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Figure 6: classical stability using Eq. I|27| l for the top panels and absolute stability using Eq. Il26p for the bottom panels. The 
3D analytical lines of Eqs. I|63l I83| l for classical stability agree with the corresponding numerical data. Above the zero-horizontal 
axes in the top panels, the Q-balls are classically unstable. Similarly, Q-balls above the horizontal axis, Eq = mQ, are absolutely 
unstable. The one dimensional Q-balls are always classically stable. The ID slopes Eq/mQ have different behaviours depending 
on DVP and NDVP unlike the other dimensional cases. 



Both DVP and NDVP analytical values uja in Table: [3] agree well with the numerical ones 
in Table: HI Generally speaking, the higher dimensional Q-balls are more stable classically as 
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well as quantum mechanically. Moreover, thin wall Q'balls are always classically stable as 
demonstrated in Eq. flU^ . but the classical stability of thick wall Q-balls is model- and D- 
dependent as in Eq. ( l83l) . The one- and two- dimensional Q-balls have a much richer structure 
than the thin and thick wall Q-balls. It is a challenging task to understand their intermediate 
profiles [7i] . 
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Table 3: Virial relations: LOa in terms of space dimension D and ratio S/IA, see Eq. 1134 l l 
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Table 4: The critical values for classical stability, absolute stability and stability against fission in DVP and NDVP using 
Eqs. ([26] |27] [28ll and Eq. The critical values are defined by ^| = = ^ (■^)| = 0' ^q/QUc. = 

and ^ = 0. The numerical values of oJa coincide with the analytic ones in Table: |3] We have confirmed numerically that 

Legendre relations Fig. [7] shows the Legendre relations: v. u, — v. Q, and v. 
iu;^ which can be used to check Eq. ( fTBj) . We have also checked the validity of the Legendre 
transformations in Eqs. fll0til2p . Since the numerical results match our analytical ones, these 
results strengthen the validity of our analytic arguments. 

5 Conclusion 

We have numerically and analytically explored the stationary properties of a single Q-ball for 
arbitrary spatial dimension D. With the time-dependent non-linear solutions in the system, 
the virial theorem induces the characteristic slopes Eq. ( l33l) . and gives the approximate critical 
values for Ua in Eq. without requiring a knowledge of the detailed profiles and potential 
forms. By linearising the Q-ball Eq. (fT9l) or rescaling in S^, we have been able to consider the 
two limiting cases called the thin and thick wall Q-balls. The step-like ansatz of Eq. (1361) can 
describe thin wall Q-balls in the extreme limit u = u^, whereas the modified ansatz Eq. (HOl) is 
applicable to ctq — cr+ which leads to wider range of parameter space u and of course includes 
the previous limit. On the other hand, the limit u ^ u+ is used to describe thick wall Q-balls 
in both the Gaussian ansatz Eq. flU7|) and our modified ansatz for the thick wall case. 



29 



DM 
0.92 
OM 



dEJdQ=tii 

.'Vumerkal d&i ill 2D + 
Numerical data in 3D 
Numerical dam in 4D 

0.92 0.94 0.9c OM 



I(l> 



10' 



Itf 



DVP 










-dS,/d<j^ - - - - 
Numerical data in 2D + 
Numerica! daai in 3D 
Numeric^ data in4D 



llf 



11/ 



dGJdl=mof 
Numerieal dukt in 2D 
.'Numerical data in 3D 
Numerieal data in 4D 



0.40 

in«f 



0.45 



IM 
OM 
0.96 
0.94 
0.92 
0.90 



dt^dQ=(0 

Numerical data iti 2d + 
Numerical data in 3D 
Numerical dam in 4D 

0.90 0.92 0.94 0.96 OM 



10> I 



10' 



^lo' 



10' 



llf 



lit 



Numencal data in 2D 
Numerical data in 3D 
Numerical dam in 4D 

lOr m' 



0.46 
0.44 



11/ 



clGJdl=mof 
Numerical dam in 2D 
.Numerical data in 3D 
Nutnencal data in 4D 



0.36 0.3S 0.40 



0.42 0.44 



0.46 0.4S 0.50 



Figure 7: The Legendre relations of Eq. JTell: = uj (left), = Q (middle), and ^ = ^uj^ (right). Note the excellent 

agreement between the analytical dotted lines and the numerical data (dots). 



The thin wall approximation is valid ioi D > 2. Since the step-like ansatz in the thin wall 
approximation does not have surface effects, the characteristic slope is simply Eq/ujQ = 1, 
Eq. ( |38l) . With the modified ansatz including surface effects, the classical stability for thin wall 
Q-balls does not depend on D in Eq. (l63l) . but the absolute stability condition Eq. ( l65l) does. 
Throughout the analysis, we have assumed Eqs. ( HTl H6l) . and imposed Eq. ( l50l) explicitly, which 
differ from the analysis in [SU]. Without these approximations, our calculations in particular 
Eqs. (13^ [31]) and Eq. fl3B|) become inconsistent. The mechanical analogies and the numerical 
results naturally explain and validate our underlying assumptions: the core sizes of the Q-balls 
are much smaller than their corresponding thickness as seen in the middle two panels of Fig. [3l 
and the surface tension depends weakly on ui as seen in Table: |2l With these assumptions, 
thin wall Q-balls for u < Ua are absolutely stable. Moreover, the characteristic slopes coincide 
with those derived using the virial theorem. This follows from our analysis of the relative 
contributions between the potential and surface energies. The slopes have two types in either 
non-degenerate vacua potentials (NDVPs) or degenerate vacua potentials (DVPs): NDVPs 
have a large energy from the charge, hence the surface energy is less effective than the potential 
energy. They support the existence of Q-matter in the extreme limit, uj = uj_. DVPs, however, 
have negligible energy from the charge compared to surface and potential energies, thus the 
surface energy is well virialised with the potential energy. As seen in the left-bottom panel of 
Fig. m the configurations of energy density have peaks within the shells, which contribute to 
the surface energy. It would be worthwhile understanding these peaks in terms of our modified 
ansatz. Even in the extreme thin wall limit, the charge and energy of the Q-balls in NDVPs 
are not proportional to the volume, i. e. no Q-matter. 

Thick wall Q-ball solutions naturally tend to free charged and massive particle solutions 
Eq. (1251) . With the simple Gaussian ansatz we have extremised S^^ with respect to ctq and R 
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with fixed cu, while the approaches in [12] are that Eq is extremised with respect to only R. By 
extremising with respect to two degrees of freedom we are able to recover the expected results 
of Eqs. (1701 [7S1) unhke in [12]. The Gaussian ansatz, however, is valid only for = 1 because 
of Eq. ( ffTl) . and gives contradictory results for the condition for classical stabihty. In order 
to remove these drawbacks in the Gaussian ansatz, we introduced another modified ansatz 
and used the Legendre relations to simphfy the computations of S'^, Q, Eq. We obtained a 
consistent classical stability condition Eq. flH^ which depends on D and a non-linear power n of 
the polynomial potential Eq. flH^ . Not surprisingly, our numerical results suggest the modified 
ansatz is much better than the Gaussian ansatz in the bottom two panels of Fig. [5l With the 
same panels, the validity condition Eq. (IHTj) in the modified ansatz has also been confirmed 
numerically. 

In Eqs. (134^ (HSD and Table: HI the analytical and numerical results found the critical value 
Ua with an assumption. The assumption says that the higher dimensional Q-balls could be 
applicable to the thin wall approximations over a wide range values of u. Although this state- 
ment may not hold for extremely flat potentials [351 because fi can be as small as 1/ Rq (see 
Eq. PUI) ) and the energy spectrum has the following proportionality Eq oc Q^^^~^^, we believe 
that it may apply to other large types of the Q-ball potentials. In summary, the higher dimen- 
sional Q-balls can be simplified into the thin and thick wall cases, while it is more challenging 
and interesting to understand stationary properties of one- and two- dimensional Q-balls. For 
example, those Q-balls embedded in 3D space (called Q-strings and Q-walls) may exist in the 
formation of three dimensional Q-balls [URl [75] . 

The properties of non-thermal Q-balls can lead to different consequences compared to ther- 
mal ones, i. e. in the evolution of the universe. The thermal effects on Q-balls induce subsequent 
radiation and evaporation. The Affleck-Dine condensate provides a natural homogeneous con- 
densate with small quantum fluctuation, these fluctuations are then amplifled to non-linear 
objects namely Q-balls if the pressure of AD condensates is negative. The formation, dynam- 
ics, and thermalization might have phenomenological consequences in our present universe, e.g. 
gravitational waves [76] and baryon to photon ratio. 
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